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Abstract 

A contemporary physicist would be hard put to agree entirely with the author of a 1959 text¬ 
book on quantum mechanics, who wrote: "A second simple, one-dimensional system, somewhat 
divorced from reality but illustrative of the principles of the theory, is a particle in a box with finite 
walls." Interest in this prototypic system has not diminished over the years, and is not likely to 
do so in the near future, because it has now been wedded to reality and is seen as a paradigm 
for quantum wells and other nanosystems. The treatment of this topic in standard textbooks has 
changed little over the last few decades, and such changes as have come to our notice leave, in 
our opinion, much to be desired. The purpose of this article is to enable, even encourage, other 
teachers to give the finite square well the attention that it deserves; to this end, the authors pro¬ 
vide a tutorial review that is more instructive and comprehensive than the accounts presented in 
numerous textbooks and dispensed in still more numerous online resources, but does not make 
any additional demands on the mathematical abilities of the student. Given the elemental nature 
of the topic, the authors claim no originality; with nothing more than the winnowing fork in their 
hands, they can only clear the threshing floor, gather the wheat into the barn, and let others decide 
whether or no the chaff should be burnt with unquenchable fire. 


1 Introduction 

When teaching the finite square well to stu¬ 
dents, teachers would do well to follow 
Lord Chesterfield's advice to his son [1, 
p. 41]: "In truth, whatever is worth do¬ 
ing at all is worth doing well, and nothing 
can be done well without attention: I there¬ 
fore carry the necessity of attention down 
to the lowest things, even to ...". We have 
omitted the last three words of his remark 
because they should be replaced, in the 


present context, by the phrase "the last sig¬ 
nificant digit". One of our aims is to reiter¬ 
ate what was pointed out (in the 1970s) after 
the advent of hand-held electronic calcula¬ 
tors [2-4]: finding the eigenenergies of a fi¬ 
nite square well provides a springboard for 
introducing a physics student to equations 
whose solutions can only be found through 
an iterative or recursive technique. To il¬ 
lustrate the operations involved in such cal¬ 
culations, we use a commercial spreadsheet 
(Microsoft® Excel®). The widespread avail- 
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ability of spreadsheet software has made 
it possible to go beyond the graphical ap¬ 
proaches which were developed in the days 
when computing facilities were primitive 
and scarce. We also show, by considering 
a specific case (a square well whose upper¬ 
most energy level is close to the top of the 
well), that occasionally the desired result 
may be obtained more easily if one is pre¬ 
pared to supplement the resources of Excel 
with some cerebration. 

We conclude the preamble by citing the 
source of the remark quoted in the Abstract 
[5, p. 41] and by recalling a recent article 
where the utilitarian relevance of the prob¬ 
lem is acknowledged [6]. 

2 Notation and Basic Relations 

Our purpose will be best served by adher¬ 
ing to the notation used in a book on mod¬ 
ern physics which has gone through three 
editions and several reprints; the symbols 
U and L stand in this book for the height 
and width of the well, respectively; m and 
E (0 < E < U) denote the mass and energy 
of the particle, respectively, and the subse¬ 
quent treatment requires the introduction 
of two more symbols [7, pp. 209-212]: 

k = (2mE/h 2 ) 1/2 , (1) 

a= [2m(U - E)/h 2 ] 1/2 . (2) 

The allowed energy levels (E^) of a 
particle confined to an infinite square well 
(of width L ) are given by the expression 

E ^° ] = = n2 ^°° 1, Say ’ (3) 

where E {°°^ is the ground state energy and n 
denotes the quantum number (n— 1,2,...). 
For the finite square well, the wavefunction 
penetrates the classically forbidden region. 


and the so-called penetration depth (or de¬ 
cay length) is defined as 


5 = 


1 

a 


h 

[2m(U — E)] 1/2 ' 


(4) 


Solution of the time-independent 
Schrodinger equation for the finite well 
leads to the conclusion that allowed values 
of E are those which satisfy the following 
transcendental equation [7, pp. 211-2]: 

(a/k) cos(kL) — sm(kL) a 

(a/k) s'm(kL) + cos(kL) k 


One sees from Eqs. 1 and 2 that both a and 
k depend on the value of E. It is convenient 
to denote the ratio a/k by T, to note that 


« lu 


~k~\E~ ’ 

(6) 

and go on to rearrange Eq. 5 as 


2 T 

tan (kL) = 2 , 

(7) 


so that one may avail oneself of the trigono¬ 
metric identity 


tan A 


2 tan(A/2) 

1 — tan 2 (A/2) ’ 


( 8 ) 


and express Eq. 7 as 

2tan£ 2T 

1 - tan 2 £ 1 - T 2 ’ U 

where 

«= f = 1^' (10) 

Now, Eq. 9 is simply a second-order 
equation in tan £ with solutions of even and 
odd parity: 


tan £ = T (Ha) 

tan £ = —1/T (11b) 


The two solution cannot coincide, for that 
would imply tan 2 £ = —1. The energy is 
now given by the relation 

= = <12 > 
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Multiplication of both sides of Eq. 5 
with the denominator on the left-hand side, 
transposition of the right-hand side of the 
resulting equation to the opposite side, and 
multiplication by k 2 leads to the result 
shown below: 


(a 2 — k 2 ) sin (kL) + 2akcos(kL) = 0. (13) 
Let us introduce the symbol 

R = VW + a 2 ) = 2AA, (14) 

h 

and use it for defining ?/ through the rela¬ 
tions 


k 

sin rj = — = 


kh 


and 


a 

COST] — — — 


R y/(2mU) ’ 

ah 


R y/(2mU)' 


(15a) 

(15b) 


which has been available to the English- 
reading public since the publication (in 
1958) of the first edition of the translated 
version of the splendid textbook on quan¬ 
tum mechanics by Landau and Lifshitz [ 8 ]. 
It is most convenient for our purpose to re¬ 
fer to the abridged version of their book [9], 
and we point out that, if our symbols L and 
U are replaced by a and U 0 , our Eq. 20 be¬ 
comes identical with their Eq. 1 on p. 80. 

We will continue to use our own nota¬ 
tion for describing the approach followed 
by Landau and Lifshitz (L&L, for short). 
They stated that their Eq. 1 can be recast in a 
more convenient form if one introduces two 
new variables, one being £ and the other 

7 = (H/LW(2/mU). ( 21 ) 

When n is even, the resulting equation is 

cos£ = ±7£, (22) 


When one recalls the identities 

2a /c 

sin(2r/) = 2 sin q cos q = ——, (16a) 

R- 

q/2 _ 2 

cos ( 277 ) = cos 2 r) — sin 2 rj = ———, (16b) 

R 1 

one sees without further ado that Eq. 13, af¬ 
ter it is divided by R 2 , can be recast as 

sin (kL + 2 rj) = 0, (17) 

or, after replacing kL by 2£, as 


and the roots for which tan £ > 0 must be taken. 
When n is odd, we have 

sin£ = ±7£, (23) 

and the roots for which tan £ < 0 must be taken. 

We will now introduce a parameter 
(called the power or strength of the well) that 
has been used by many authors [10-16]; we 
will denote it by the symbol P, and it is 
enough to point out that P = I/ 7 , and to 
note, for later use, that 


sin [2 (£ + 77 )] = 0 . (18) 

For reasons that will become transparent 
shortly, we observe that Eq. 18 implies that 

2£ = (n + 1)tt - 2r/ (19) 



(0 < £ < P ). (24) 


Finally, putting E = 0 in Eq. 4, and denot¬ 
ing the resulting expression by A, we write 


where n = 0,1, 2, • • •. Replacing 2£ by kL A 

and rj by sin -1 [kh/y/ ( 2 mU)\, we arrive at the 

relation and note that 

A 

kL — (n + l) 7 r — 2sm~ 1 [kh/y/(2mU)], (20) L 


h 

V(2 mUy 

1 

2 ' 


(25) 

(26) 
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3 A Retrospect 

The first edition (in Russian) of L&L's 
book on non-relativistic quantum mechan¬ 
ics was published in 1947; its English trans¬ 
lation appeared in 1958 [8]. GoTdman, 
Krivchenkov, Kogan and Galitskii (here¬ 
after GKKG) used L&L's approach but 
added a graphical solution in their well- 
known book, whose English translation 
was published in 1960 [11]. Among the au¬ 
thors who wrote directly in English, Pitka- 
nen [10] appears to have been the first 
to propose some graphical representations 
based on Eqs. 22 and 23, and one of these 
representations is equivalent to that dis¬ 
played in Fig. 20 of Ref. [11, p. 60]; it is 
not inconceivable that, en route, he also de¬ 
duced our Eq. 20, but found it unnecessary 
to bring it to the readers' attention because 
in those days teachers, more occupied with 
finding user-friendly graphical procedures, 
spent a great deal of ingenuity on devising 
plots which used simple continuous curves 
(such as sinusoids) and straight lines. In 
an age when a contemporary student can 
plot, with a few keystrokes and a couple of 
mouse clicks, half a dozen curves, each with 
a different colour, the relevant question to 
ask is not which graph is easier to gener¬ 
ate but which is more instructive than other 
claimants to a student's attention. 

Cantrell [17], who referred to Landau 
and Lifshitz [8], concluded his contribution 
with the following comment: 

The only pitfall of the method is in remember¬ 
ing which diagram (sine or cosine) is appropriate 
for a given situation, and where the "forbidden" 
regions are on the diagram. While a complete 
mnemonic would be too lengthy to be useful, it 
is possible to give a simple rule; the sine, which is 
an odd function, is associated with the odd-parity 
states in the symmetrical well..., and with nega¬ 
tive values of tan ka [tan(fcL), in our notation], for 
the purposes of this construction. 


Now, spoon-feeding is usually defined as 
'providing someone with so much help or 
information that they do not need to think 
for themselves' [18]. Students who are 
unable to figure out, after thinking hard 
(and, if necessary, long) the restrictions on 
the sign of tan £ should be encouraged to 
change to a field of study better suited to 
their particular talents. A key reference 
in our bibliography [13] has a space filler 
in which Sylvia Plath is quoted as saying, 
"The day I went into physics class was 
death"; one only hopes that the depres¬ 
sion which tormented the prodigiously tal¬ 
ented Plath has not been traced by a biogra¬ 
pher to tiresome mnemonics composed by 
a physics teacher desperate to narrow the 
gulf between "The Two Cultures". 

An equation of the same form as Eq. 13 
was derived by Reed [16], and subse¬ 
quently shown by Sprung, Wu and Mar- 
torell (SW&M) to be equivalent to the com¬ 
pact form displayed in Eq. 18; in their 
first publication on the topic [14], SW&M 
claimed (justifiably) to have derived no 
more than a "new and very accurate ana¬ 
lytic solution ... for the energy levels of the 
finite square well potential, in the form of a 
rapidly convergent series in inverse powers of 
the strength" (emphasis added). The discov¬ 
ery of Eq. 18 was attributed by Aronstein 
and Stroud [19] to SW&M [14], who began 
the sequel to their first paper with a claim 
for the derivation of the simple formula it¬ 
self [15]: 

Some time ago we provided a simple analyti¬ 
cal formula for the bound state energies of the fi¬ 
nite square well potential. We also showed how 
our solution is related to the graphical construc¬ 
tion of Pitkanen. We believe that these meth¬ 
ods are superior to those presented in most text¬ 
books, and should help in presenting this canon¬ 
ical example in introductory quantum mechanics 
classes. 
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We conclude this section by pointing out 
that explicit analytical expressions for the 
energies of the finite square well are also 
available [20-23], but will not be consid¬ 
ered here, because their derivation and fur¬ 
ther manipulation requires techniques not 
found in the mathematical repertoire of a 
physics undergraduate; we will, however, 
quote and use an approximation that is 
valid for the case of small 7 (see below). 


4 Calculation of Eigenenergies 

Since we wish to comment on the works of 
other authors, we will use two conventions 
for actual calculations, and measure ener¬ 
gies in eV and distances in nm, or work 
with dimensionless quantities; in the latter 
case, energies will be measured in terms of 
£ = h 2 /mL 2 and lengths in terms of L. 


4.1 First Choice: Eq. 7 

If one's sole aim is to find the eigenenergies, 
one need go no further than Eq. 7, since it 
is a single equation in a single unknown 
E; in a graphical representation, one would 
look for the intersections of the curve y = 
tan( A'L), which shows multiple discontinu¬ 
ities, with the curve y = 2r/(l — T 2 ), which 
has only one singularity (at E = U /2). 

Let us rewrite Eq. 7 in terms of the sym¬ 
bols representing the properties of the par¬ 
ticle and the well. 


2y/(U-E)E 

2E-U 


tan 



(27) 


and consider a particular case that was an¬ 
alyzed recently by Lindberg [24]: a well of 
height U = 25 eV and L = 0.5 nm. Lind¬ 
berg prepared a three-column table and put 
values of energy E in steps of 0.1 eV from 
0.1 to 25 eV (in the first column) and the 
corresponding values of the two sides of 


Eq. (27) in the other two columns. He 
plotted the values of the second and third 
columns against E, looked for intersections 
of the two curves, and adjusted the energy 
at each intersection "to get something close 
to a perfect fit, i.e. LHS = RHS". 

The procedure followed by Lindberg is 
simple and robust, and adequate if the 
teacher's aim is to obtain results correct to 
a few decimal places. A minor blemish is 
that the plot, which exhibits five crossings, 
does not separate levels belonging to differ¬ 
ent parities, since the first three crossings 
occur when the left-hand side is negative 
(0 < E < U/2), and the last two when the 
left-hand side is positive (for U/2 < E < U). 

We believe that the pedagogical value of 
the exercise can be substantially enhanced 
by informing the students that other, more 
efficient ways of finding precise values of 
the roots are available. We will describe 
these methods below before considering 
the advantages of using other equations for 
finding the eigenenergies. 

4.2 Some Root-Finding Procedures 

The Analysis Tool Pa k of Excel contains two 
tools, namely Goal Seek and Solver [25], 
which can find (in most cases) the roots of 
the equation f(x) = g(x) within a stated 
precision; the user must supply an initial 
guess for the root, say x 0 , calculate the value 
of the difference 8(x 0 ) = f(x 0 ) — g(x 0 ), and 
let the tool search for the value of the argu¬ 
ment for which 8 attains a value sufficiently 
close to zero (that is, less than a value set 
by the user). Goal Seek and Solver are su¬ 
perb tools indeed, but it seems advisable to 
inform a novice, who is likely to view each 
of these as a magic wand, that they are not 
foolproof (see below) and other methods, 
not making explicit use of Analysis ToolPak, 
are available. 

In the mid-1970s, when personal com¬ 
puters were not yet commonplace, the 
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availability of hand-held electronic calcula¬ 
tors (also called electronic slide rules) was 
seen by some teachers as an opportunity 
to introduce their students to root-finding 
procedures. We will refer to three brief 
and lucid contributions [2-4], which fo¬ 
cussed primarily on the finite square well, 
and will take it for granted that those who 
need some instruction would consult these 
sources, and other references cited therein. 
Among the many methods that could be 
brought to bear on our problem, only two 
will be considered here, namely iteration 
and bisection. 

The basic idea behind the bisection 
scheme has been explained in the context 
of the finite square well in a recent book 
on computational physics [26, pp. 134-135]; 
the Internet is a rich source of information 
as to how this method may be implemented 
by using Excel, and it will be enough to cite 
one resource here [27], which has tested by 
the authors and found to be satisfactory. 

As it stands, Eq. 7 is not of the form 
x = f(x), and thus not well-suited to an it¬ 
erative treatment, but it can be easily trans¬ 
formed, as shown above, in two different 
ways into a pair of equations with the de¬ 
sired form; details concerning an iterative 
extraction of the roots of the transformed 
equations are given below. 

The five allowed energy levels of the 
square well under consideration are dis¬ 
played in Table 1. The first four roots of 
Eq. 27 found by using bisection. Goal Seek 
and Solver concurred with one another. It is 
crucial to note that Solver failed to find E 5 , 
regardless of how close the initial guess (in 
eV), to be denoted as El, was to the cor¬ 
rect value, but Goal Seek and bisection re¬ 
turned the correct answer, so long as the in¬ 
put from the user was chosen judiciously; 
as to what judicious choice means in this con¬ 
text, there is no substitute for personal ex¬ 
perience, but a few remarks may not be 
amiss here. Though we will speak of only 


Goal Seek, the following remarks are equally 
applicable to Solver. 

Table 1: Energy levels of a finite square well 
(U = 25 eV, L = 0.5 nm). The calculations 
used NIST data [28], which led to the following 
values: he = 1.97327 x 10 2 eV-nm and me 2 = 
5.10999 x 10 5 eV. 


n 

EjeV 

1 

1.12294 

2 

4.46186 

3 

9.90751 

4 

17.16578 

5 

24.78411 


Let us consider the third eigenvalue, and 
note that if one takes El = 9.3, which is 
much closer to E :i than to E 2 , Goal Seek con¬ 
verges to the value of E 2 , furthermore, for 
many choices of El, for example 9.4 or 12.4, 
Goal Seek fails to find any of the eigneval- 
ues. An alert student should be able to 
explain these observations (and others of 
a similar nature) by examining a graph in 
which y = 2a /(U — E)E/ (2 E — U) and y = 
tan y/(2mEL 2 /h 2 ) are plotted against E. 

4.3 Separating States of Even and Odd 
Parity 

As already stated, to anyone who has access 
to modern computing facilities, the differ¬ 
ence between plotting two curves or three 
is not an important consideration. It is 
worthwhile, therefore, to employ a differ¬ 
ent graphical procedure for illustrating the 
calculation. Figure 1 is a plot based on 
Eqs. 11a and lib; y = tan£, y — T and 
y = —1/T have been plotted against £. One 
sees from the figure that three solutions of 
even parity and two of odd parity are per¬ 
mitted, and that the even/odd parity solu¬ 
tions correspond to the intersections in the 
upper/lower half of the plot. 
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0 1 2 3 4 5 6 

5 


Figure 1: Graphical representation of Eqs. 11a and lib (U = 25 eV, L = 0.5 nm). The intersections 
of y = tan £ and y = T, labelled as 1, 3, and 5, belong to even parity states; the intersections of 
y = tan£ and y = — 1/T, labelled as 2 and 4, pertain to odd parity states. The vertical arrow 
indicates the abscissa (£ max = P = 6.4039 • • •) at which y = T crosses the horizontal axis. The 
dotted vertical line is the asymptote to y = tan £ at £ = 7t/2. 


The pivotal role played by the parame¬ 
ter 7 (or its reciprocal P ) becomes appar¬ 
ent when one considers two limiting cases: 
P — y oo and P —* 0. A student ought to be 
able to deduce, merely by examining the lo¬ 
cations of the intersections, that in the for¬ 
mer case, one recovers Eq. 3, and that if 
P < 7t/2, only one energy level of even par¬ 
ity will be allowed. We suggest that the task 
of supplying the missing steps in some re¬ 
marks made by L&L and GKKG may also 
be used as student exercises. In treating the 
case 7 <§; 1, L&L argue as follows: "In par¬ 
ticular, for a shallow well in which U 0 ££> 
h 2 /ma 2 , we have 7 C 1 and equation (3) 
[our Eq. 23] has no root. Equation (2) [our 
Eq. 22] has one root (with the upper sign on 
the right-hand side), £ = I/7 — I/27 3 . Thus 
the well contains only one energy level,..., 
which is near the top of the well." GKKG 
also made this remark (without an explana¬ 
tion), and went on to state that the number 


of levels for arbitrary values of U and L will 
be equal to N, where N follows from the re¬ 
lation 

2 P 

N> — >N- 1. 

7r 


In modern notation the above result is 
stated as 


N = 1 + 


P 

yj2 


(28) 


where [xj is the floor function defined here 
as the largest integer smaller than x. 

We turn next to the problem of finding 
the abscissa for each intersection, and will 
consider four options: Goal Seek, Solver, bi¬ 
section and iteration. Phillips and Murphy 
have discussed the solution of the equation 
x = tan(/3x) by iteration [2], and their dis¬ 
cussion can be easily adapted to Eqs. 11a 
and lib. It will be enough to state here that 
Ei and E 3 were determined by transform¬ 
ing Eq. 11a as £ = tan -1 (T), and Eq. lib was 
rephrased as £ = tan _1 (—1/T) for deter- 
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mining E 2 and E 4 ; for the uppermost level, 
T was considered to be the independent 
variable, and Eq. 11a treated as T = /(T). 
All four methods gave results which were 
concordant, except that Solver failed to find 
the value of E 5 (which is exceptional in the 
sense that it lies close to the top of the well). 

The analysis presented above was devel¬ 
oped by the authors as an alternative to the 
use of Eqs. 22 and 23, but other teachers 
may prefer to use the latter pair. Appendix I 
deals with the application of the iterative 
approach to these equations. 


5 Finite and Infinite Wells 


Since the wavefunction of a particle subject 
to a finite square well potential extends be¬ 
yond the boundaries [7, pp. 209-212], de¬ 
caying exponentially at a rate 5, Garrett 
[12] proposed that the approximate energy 
eigenvalues (to be denoted by £ n ) of the fi¬ 
nite square well can be approximated by 
adapting Eq. 3, valid for the energies of an 
infinite square well, as follows: 


£ = 


n 2 n 2 h 2 

2m(L 4- 2rL') 2 ' 


(29) 


Since 5 n is an energy-dependent quantity, 
Eq. 29 is, in fact, an implicit relation for £ n 
that must be solved numerically for a given 
value of n by using the relations 



h 



U - £ { n~ X) ] 


with k = 

£(k) 


1 ,2,..., and 
n 2 n 2 h 2 


2m 


L + 2 5 { n k) 


2 ’ 


(30) 


( 31 ) 


where £n ! and d (k) each denote the A -th it¬ 
erate. All that is needed now is an appro- 


§ These words are form T. S. Elliot's Burnt Norton. 


priate guess for £n\ and Garrett suggested 
that one should take = E^°\ 

The authors of Ref. [7], hereafter re¬ 
ferred to as SM&M, illustrated Garrett's ap¬ 
proach by calculating £\ k] for a well with 
U = 100 eV and L = 0.2 nm. To start 
the iteration, they set £ i \ )) = 0 (not = 
E^°\ as proposed by Garrett!), and ob¬ 
tained thereby 

fc 

6™ = = 0.0195 nm, (32) 

V2 mU 

and 


£[ 1] = 


7T 


■K 2 


2m 


L + 25\ 


(O' 


n2 


= 6.58 eV. 


(33) 


They went on to calculate 




( 2 ) 


h 



U-£ 


(i) 

i 


2 


0.0202 nm, (34) 


inserted it into Eq. 31 to find £ (2> = 

6.53 eV, and stated: "The iterative process 
is repeated until the desired accuracy is 
achieved. Another iteration gives the same 
result to the accuracy reported. This is in 
excellent agreement with the exact value, 
about 6.52 eV for this case." 

The above calculation and the ensuing 
conclusions are so farcical as to be fit only 
for a homework exercise designed to test 
a student's ability for plugging numerical 
values into a formula and calculating the 
answer. A student who is given this task 
would be able, if not incapacitated by "tu¬ 
mid apathy and no concentration",^ to re¬ 
produce the values of <5^, £[ v) and given 
above (Eqs. 32-34), and would probably 
mutter "So far so good"; if the student has 
set up the spreadsheet formulas correctly, 
(s)he would find that £[ 2) = 6.5071 eV, and 
£{ 3) = 6.5080 eV. The discrepancy between 
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these values and the results of SM&M (ac¬ 
cording to whom £j 2) = 6.53 eV = S[ 3 \ 
and E x = 6.52 eV) might prompt a stu¬ 
dent (or a teacher) to calculate E 1 for her¬ 
self/himself, and such a person would find 
that Ei = 6.557 eV. If one repeats the cal¬ 


culation and takes £-[ 0) = E^°\ one would 
find £ j 2) = 6.508 eV, which shows that con¬ 
vergence is rather rapid in this particular 
case, and that, whatever the initial guess, 
one ends up with a result that is close to 
6.508 eV. The results of our own calculation 
have been reproduced in Figure 2. 


n 

Finite Square Well 

Infinite Square Well 

i 

£' 1 =6.55679 eV 


i? 1 [ °° 1 =9. 40075 eV 


i 


5^ 

<G (i) 


8 1 (,) 

<G (i) 

l 

0 

0.01952 

6.58093 

9.40075 

0.02051 

6.47350 

2 

6.58093 

0.02019 

6.50714 

6.47350 

0.02018 

6.50839 

3 

6.50714 

0.02019 

6.50800 

6.50839 

0.02019 

6.50799 

4 

6.50800 

0.02019 

6.50799 

6.50799 

0.02019 

6.50799 

5 

6.50799 

0.02019 

6.50799 

6.50799 

0.02019 

6.50799 

6 

6.50799 

0.02019 

6.50799 

6.50799 

0.02019 

6.50799 

7 

6.50799 

0.02019 

6.50799 

6.50799 

0.02019 

6.50799 

8 

6.50799 

0.02019 

6.50799 

6.50799 

0.02019 

6.50799 


n 

Finite Square Well 

Infinite Square Well 

2 

is 2 =12.2471s 

E 2 M = 19.73921s 


i 


§2 ( ° 

<G (0 


§2 (<) 

£P 

1 

0 

0.12500 

12.63309 

19.73921 

0.20194 

10.01539 

2 

12.63309 

0.16068 

11.30553 

10.01539 

0.15081 

11.65101 

3 

11.30553 

0.15544 

11.48699 

11.65101 

0.15675 

11.44107 

4 

11.48699 

0.15612 

11.46298 

11.44107 

0.15595 

11.46908 

5 

11.46298 

0.15603 

11.46617 

11.46908 

0.15606 

11.46536 

6 

11.46617 

0.15605 

11.46575 

11.46536 

0.15604 

11.46586 

7 

11.46575 

0.15604 

11.46580 

11.46586 

0.15604 

11.46579 

8 

11.46580 

0.15604 

11.46580 

11.46579 

0.15604 

11.46580 


Figure 2: Excel output showing the application 
of Garrett's approximation to a square well with 
U = 100 eV and L = 0.2 nm. The first three 
columns show the results obtained by using 
s[ 0) = 0; the last three, by using £-[ 0) = 


Figure 3: Excel output showing the application 
of Garrett's approximation to a square well with 
U = 32e, where e = h 2 / (mL 2 ). The first three 
columns show the results obtained by using 
d) 01 =0; the last three, by using £< 2 (>) = 


Garrett (who used the symbol a for our 
L) considered a well of depth 32 e, put 

£-[ 0) = E^°\ and he too stopped at the sec¬ 
ond iteration, not only for calculating the 
energy of the first level, where convergence 
is indeed rapid, but also for calculating £ 2 . 
He found the exact energy to be E 2 = 12.0 e, 
took £f = Til, 00 = 2 tt 2 c, and arrived at 
the result £ 2 2) = 11.65 e. The good agree- 
ment between £ 2 and E 2 is illusory. One 
can verify, in the first place, that the exact 
value reported in Ref. [13], E 2 = 12.2471 e, 
is correct, and (secondly) that, if the iter¬ 
ation is continued, convergence is reached 
after seven steps (see Figure 3) and yields 
£ 2 " ) = 11.4658£ (m > 7). It is curious that 


the authors of Ref. [13], who were interested 
in comparing the relative error in their first- 
order approximation with that in Garrett's 
approximation, did not press the latter ap¬ 
proximation to its logical end. 

The well considered by Garrett has only 
three bound states. He stated: "Clearly, this 
method cannot be used for the n = 3 state 
since E 0 > Go [G'°° ] > U, in our notation], 
... The value of £^ in [his] Table 1 was ob¬ 
tained by taking a linear extrapolation of 5[ 
and 5' 2 to get <5 3 = 0.171a". Garrett's sym¬ 
bols 5[ and S 2 correspond to our ()[' j and 
8^, respectively, and his <5 3 , which corre¬ 
sponds to our , has no superscript be¬ 
cause he did not continue the iteration. This 
stratagem led him to the result £^ ] = 24.7 e, 
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in excellent agreement with the exact value 
E 3 = 25.9 eV. If Garrett had proceeded to 
the next step of iteration, using as the 
input, he would have found = 19.IT e:! 
Since the iteration under consideration al¬ 
ways converges to the same result for any 
choice of £n ] that is lower than U (see Fig¬ 
ure 3), the only legitimate result is that 
found after convergence is reached ; one must 
conclude that Garrett's approximation does 
not provide the accuracy he claimed. 

Although Garrett's scheme does not pro¬ 
vide useful results, the idea of replacing a 
finite square well by a wider infinite well 
can be placed on a sound basis. The reader 
will have noticed that Garrett's own pro¬ 
posal amounts to replacing a finite square 
well (A allowed energy levels) with N in¬ 
finite wells, each of a different width. An 
obvious amendment is to replace the finite 
well with a single infinite well of a unique 
width L; to this end, we remind ourselves 
of Eqs. 4, 25 and 26, and define 


and 


L — L + 2A — T(1 + y), 
- n 2 ir 2 h 2 

' J n 

2 mL 2 


(35) 

(36) 


One expects that the approximation would 
work well when E n <C IJ, and that it would 
gradually worsen as n becomes larger, and 
this is indeed found to be the case [13]. 

Barker and co-authors [13] noted that a 
Taylor expansion of Eq. 22 about mr/2. 



sin(n7r/2) 

1 ! 



sin(wr/2) 

3! 



(37) 


where n is an odd integer, could be used 
to develop approximations for Eq. 22, and 
likewise for Eq. 23. If one retains only the 
first-order term, one arrives at the result 


in — (utt/2)[P/ (P + 1)], (38) 


which leads, by virtue of Eq. 12, to Eq. 36. 

Paul and Nkemzi [21] showed that the 
explicit expression for E n can be expanded 
as a series in powers of 7 ; when their result 
is corrected for an algebraic error, identified 
by Aronstein and Stroud [22], the expres¬ 
sion for the energy comes out to be 


E n (l) 


n 2 K 2 7 T 2 
2 mL 2 


+ 0 ( 7 3 ). 


(39) 


An altogether different motivation 
(the validity of Ehrenfest's theorem) led 
Rokhsar [29] to the relation L = L{ 1 + 7 ). 

The results displayed in Figs. 2 and 3 do 
not corroborate the conclusion, reached in a 
recent publication [30], that Garrett's itera¬ 
tion, if carried ad infinitum, leads to the re¬ 
sult stated in Eq. 36. This can only happen 
if S ( l f c> becomes independent of k and equals 
A, which requires fulfillment of the absurd 
condition £jf° ! = 0 . 


6 Concluding Remarks 

The reader who is totally satisfied with the 
sterling treatment of the finite square well 
problem by L&L, and has also kept abreast 
of the literature, deserves our apologies, for 
nothing is more tedious than reading an 
account of the errors one has not commit¬ 
ted. Our motivation for writing this article 
has been to present an alternative graphi¬ 
cal representation; to stress that the exer¬ 
cise of finding the eigenenergies of the fi¬ 
nite well scrutinized by us will serve as 
an antidote to the tyro who may otherwise 
be tempted into thinking of Goal Seek and 
Solver as magic wands; and to suggest that 
those who teach (or use) iteration should 
obey the commandment "Thou shalt not 
terminate (kill?) an iteration prematurely". 
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Appendix I: Iterative Solution of the Equations of Landau and Lifshitz 

Since all necessary information is readily available [3], we will only state that, for the 
problem of a well with U = 25 eV and L — 0.5 nm, the five roots (<f/„ k =1-5) of Eqs. 22 
and 23 were determined by transforming these relations into those listed below: 

f = /i(0 = c°s -1 (7£), (22:1) 

£ = / 3 (f) = cos -1 (7£) + 7r, (22:3) £ = / 2 (C) = sm _1 (70 + n, (23:2) 

£ = /s(0 = cos(0/7; (22:5) £ = / 4 (f) = sin -1 ^) + 2vr. (23:4) 

The output of the Excel worksheet prepared by the authors is displayed in Figure 4. 


i 

p ('-1) 

SI 

fr 


S2 



p o-i) 

S3 



p 

S 4 



p 

S5 


1 

5 

0.67490 


5 

2.24570 


5 

3.81650 


5 

5.38729 


5 

1.81656 

2 

0.67490 

1.46521 


2.24570 

2.78330 


3.81650 

4.07393 


5.38729 

5.28361 


1.81656 

-1.55807 

3 

1.46521 

1.33995 


2.78330 

2.69197 


4.07393 

4.02288 


5.28361 

5.31289 


-1.55807 

0.08149 

4 

1.33995 

1.36000 


2.69197 

2.70775 


4.02288 

4.03317 


5.31289 

5.30475 


0.08149 

6.38271 

5 

1.36000 

1.35680 


2.70775 

2.70503 


4.03317 

4.03110 


5.30475 

5.30702 


6.38271 

6.37227 

6 

1.35680 

1.35731 


2.70503 

2.70550 


4.03110 

4.03152 


5.30702 

5.30639 


6.37227 

6.37856 

7 

1.35731 

1.35723 


2.70550 

2.70542 


4.03152 

4.03143 


5.30639 

5.30656 


6.37856 

6.37485 

8 

1.35723 

1.35724 


2.70542 

2.70543 


4.03143 

4.03145 


5.30656 

5.30651 


6.37485 

6.37707 

9 

1.35724 

1.35724 


2.70543 

2.70543 


4.03145 

4.03145 


5.30651 

5.30653 


6.37707 

6.37576 

10 

1.35724 

1.35724 


2.70543 

2.70543 


4.03145 

4.03145 


5.30653 

5.30652 


6.37576 

6.37654 

11 

1.35724 

1.35724 


2.70543 

2.70543 


4.03145 

4.03145 


5.30652 

5.30653 


6.37654 

6.37607 

12 

1.35724 

1.35724 


2.70543 

2.70543 


4.03145 

4.03145 


5.30653 

5.30653 


6.37607 

6.37635 

13 

1.35724 

1.35724 


2.70543 

2.70543 


4.03145 

4.03145 


5.30653 

5.30653 


6.37635 

6.37619 

14 













6.37619 

6.37628 

15 













6.37628 

6.37623 

16 













6.37623 

6.37626 

17 













6.37626 

6.37624 

18 













6.37624 

6.37625 

19 













6.37625 

6.37625 

20 













6.37625 

6.37625 


Figure 4: Excel output showing iterative calculations of the roots of Eqs. 22 and 23. 
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